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WIETHODE DE SIMULATION POUR PREDIRE EN FONCTION DU TEMPS UNE COMPOSITION DETAILLEE 
D UN FLUIDE PRODUIT PAR UN RESERVOIR. 


(57) Methode de simulation pour pr6dire en fonction du 
temps une composition detailt§e & Q composants et/ ou 
pseudo-composants d'un tluide prodult par un reservoir. 
Elle est du type consistant t: 

- representer le reservoir sous la forme d'un reseau de 
mailles 0) dont chacune constitue un volume elementaire 
rempli de fluide, 

- realiser, de fagon connue en sol, une simulation com- 
positionnelle du fluide regroups & N composants et pseudo- 
composants (i), N etant inferleura Q, et elle est caract6risee 
en ce qu'elle consiste en outre 

e) determiner, k chaque pas le temps (m) et pour cha- 
que maille (j) la composition des phases liquide et vapeur du 
fluide d6taille de !a composition Q du fluide, ^ partir des va- 
leurs de la fraction vaporis§e (6]"") et des constantes d'6qui- 
libre (kjO du fluide regroups, 

f) ^valuer pour chaque maille, au pas de temps (m+1 ), la 

auantit6 molaire de chacun des Q composants et/ ou pseu- 
o-composants du fluide d6taill6 h partir des valeurs corres- 
pondantes au pas de temps (m), des d6bits des phases du 
fluide regroupe et des compositions des phases liquide et 
vapeur du fluide d^taill^ d^termin^es d t'^tape e), et k 

g) ^valuer, pour chaque puits de production, la compo- 


sition detalll6e du fluide produit entre des Instants t et t' cor- 
respondant ^ux pas de temps m et m+1 , ^ partir des debits 
des phases du fluide regroupe et des compositions des pha- 
ses liquide et vapeur du fluide detalll6 determin^es k I'etape 
e). 
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Methode de simulation pour predire en fonction du temps 
une composition detaillee d*un fluide produit par un reservoir 


La presente invention conceme une methode de simulation pour 
predire, en fonction du temps, la composition detaillee d'un fluide produit 
par un reservoir et plus particulierement la composition detaiUee d'un fluide 
contenu dans et produit par un gisement petrolier dans lequel sont implantes 
10 Tin ou plusieurs puits de production. 

La simulation compositionnelle d'un gisement petrolier est 
couramment utilisee pour foumir les profils previsionnels de production du 
gisement qui permettent notamment de detenniner le schema de production 
le mieux adapte a ce gisement. 
15 La simulation compositionnelle d'un gisement petrolier est mise 

en oeuvre non pas en utiHsant une description reelle du fluide mais en 
utilisant un fluide modelise par un nombre de composants plus reduit que le 
nombre de composants du fluide reel. En effet, le nombre de composants du 
fluide de gisement etant relativement grand, la modelisation avec tons les 
20 composants conduirait a des temps de calcul trop importants. Ce temps de 
calcul prohibitif a amene les specialistes a regrouper les composants purs en 
pseudo-composants, par exemple des pseudo-composants regroupant Vazote 
et le methane, des pseudo-composants regroupant les hydrocarbures en C3 
et C4, . . . etc. et a realiser la simulation de gisement sur un nombre reduit 
25 N de composants purs et de pseudo-composants. Une telle composition 
reduite est designee par Texpression "composition regroupee du fluide" ou 
"fluide regroupe". Generalement N represente 5 a 8 composants et pseudo- 
composants, ce qui est considere comme etant suffisant pour bien 
representer le comportement du fluide de gisement aux conditions de fond. 
30 Avant d'effectuer la simulation compositionnelle de gisement, on 

represente le reservoir ou gisement sous la fonne d'xm reseau de mailles 
dont chacune constitue un volume elementaire dudit reservoir. Le nombre 
de mailles peut atteindre plusieurs milliers et chaque maille presente des 
proprietes qui lui sont propres, comme la geometric, la porosite, la 
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permeabilite, etc., De plus, a certaines mailles peut correspondre au moins 
un puits d'injection ou de production qxii est implante dans ledit reservoir. 

La simulation de gisement peimet de calculer pour chaque maiUe 
un certain nombre de variables principales du fluide de gisement, lesdites 

5 variables pouvant etre la quantite (nombre de moles) de fluide, la compo- 
sition ou fraction molaire de chaque composant et pseudo-composant du 
fluide regroupe et la pression regnant dans chaque maille. Ces variables 
sont connues a I'instant t=0 (debut de I'exploitation du gisement) et sont 
ensuite calculees par la simulation de gisement a chaque instant t ou pas de 

10 temps m. Pour chaque maille, il est possible, a partir de ces variables 
principales, de calculer toute autre propriete du fluide present dans la 
maille, comme le nombre de phases, la composition de chaque phase, etc. 
Pour chaque maLQe a laquelle est associe un puits de production, il est 
possible, a partir de ces variables principales et des contraintes de 

15 production imposees, de calculer egalement le debit de production et la 
composition du fluide produit par ledit puits. 

Les proprietes thermodynamiques du fluide peuvent etre 
calculees en utilisant une des equations d'etat bien connues des specialistes 
et qui ne seront done pas rappelees. 

20 La modehsation a N composants, si elle est sufQsante poxu* 

representer le comportement du fluide de gisement aux conditions de fond, 
n'est plus appropriee pour simuler le comportement du fluide dans les 
installations pour son exploitation en surface, laquelle exploitation necessite 
de connaitre une composition plus detaillee a Q composants et/ou pseudo- 

25 composants du fluide produit par le gisement, Q etant superieur a N et par 
exemple de I'ordre de 20 a 30. 

Jusqu'ici, on obtenait cette composition detaillee a Q composants 
et/ou pseudo-composants a partir de la composition regroupee N en 
supposant que la composition des pseudo-composants restait constante dans 

30 le temps. 

Une teUe maniere de proceder est source d'erreurs car la 
composition des pseudo-composants varie dans le temps, notamment en 
fonction de la pression du gisement ou suite a I'injection d'un gaz dans le 
gisement. 
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De ce fait, la composition detaillee obtenue reste tres approxima- 
tive et ne permet pas de predire de maniere satisfaisante le comportenient du 
fluide dans les conditions d'exploitation en surface bien des annees apres le 
debut de la production. 

La presente invention a pour but de proposer une methode de 
simulation pour predire la composition detaillee du fluide produit par un 
reservoir en tenant compte des parametres du reservoir et qui soit beaucoup 
plus precise. 

L*invention a pour objet une methode du type consistant a : 

a) representer le reservoir sous la forme d'un reseau de mailles (j) 
dont chactme constitue un volume elementaire rempli de fluide, 

b) definir le fluide par une modelisation regroupee a N 
composants et pseudo-composants (i) et a determiner une equation d'etat 
decrivant le fluide dans cette modelisation regroupee, 

c) defmir egalement le fluide par une modelisation detaillee a Q 
composants et/ou pseudo-composants, Q etant superiem- a N, et a deter- 
miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa9on connue en soi, une simulation composi- 
tionneUe du fluide regroupe a N composants et pseudo-composants (i), 
ladite simulation compositionnelle permettant de calculer au moins pour 
chaque maille (j) et a des pas de temps consecutifs (m, m+1, ...) la fraction 
vaporisee (g^), les constantes d'equilibre liquide-vapeur [k^) de chaque 
composant (i), les debits d'injection ou de production (s;;) et pour chaque 
paire de mailles (j, h) les debits des phases liquide (u^J et vapeur (u^J du 

fluide a N composants et pseudo-composants, et eUe est caracterisee en ce 
qu'elle consiste en outre a: 

e) determiner, a chaque pas de temps (m) et pour chaque maille 
(j) la composition des phases hquide et vapeur du fluide detaille de la 
composition Q du fluide, a partir des valeurs de la fraction vaporisee (ef ) et 
des constantes d'equilibre (k") du fluide regroupe, 

f) evaluer pour chaque maille, au pas de temps (m+1), la quantite 
molaire de chacun des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs coirespondantes au pas de temps (m), des debits 
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des phases du fluide regroupe et des compositions des phases liquide et 
vapetu* du fluide detaiUe determinees a I'etape e), et a 

g) evaluer, pour chaque puits de production, la composition 
detaillee du fluide produit entre des instants t et f correspondant aux pas de 
5 temps m et m+1, a partir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
Tetape e). 

Selon ime autre caracteristique de la presente invention, les 

etapes e, f et g sent mises en oeuvre en meme temps que la realisation de la 
1 0 simulation compositionnelle du fluide regroupe. 

Selon une autre caracteristique de la presente invention, les 

resultats de I'etape d sont memorises puis utilises ulterieurement pour la 

mise en oeuvre des etapes e, f et g. 

La methode selon la presente invention est mise en oeuvre pour 
15 un reservoir constitue, par exemple, par un gisement petrolier. De maniere 

cormue, on realise une simulation de gisement. A cet effet, on represente le 

gisement sous la forme d'un reseau maille, certaines des mailles ou im 

groupe de mailles etant associes a un des puits de production rmplantes dans 

le gisement petrolier a exploiter. 
2 0 Pour rediiire le temps de calcul, la simidation compositiormeUe 

du gisement est effectuee sur im nombre limite N de composants et pseudo- 

composants, par exemple de 5 a 8, defmis de la maniere indiquee 

precedenunent, ces composants et pseudo-composants etant selectioimes en 

fonction de la nature du gisement. 
25 Le flmde regroupe a N composants et pseudo-composants est 

decrit par une equation d'etat qui pourrait etre par exemple celle de PENG- 

ROBIMSON adaptee a ladite composition regroupee. 

Au debut de I'exploitation du gisement, on coimait, par des 

mesures effectuees au prealable, pour chaque maille j du reseau, les 
30 variables principales a I'instant t=0 ou pas de temps m=0. Ces variables 

principales sont ; 

la quantite de fluide : ff 

la fraction molaire du composant i : z" avec 1 ^ i ^ N 
la pression : p^ 
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PoiiT passer du pas de temps m au pas de temps m+1, la 
simulation compositioimelle comprend plusieurs etapes detaillees ci-apres. 

Dans tme premiere etape et poiir le pas de temps m, un calcul 
d'equilibre liquide-vapeur (flash) est effectue sur le fluide regroupe, pour 
5 chaque maille j afm de determiner : 

la fraction vaporisee 0 J' 

la fraction molaire de chaque composant i dans la phase huile (x" avec i = 1...n) et 
la fraction molaire dudit composant i dans la phase gaz (y" avec i = 1. . .n). 

Ces diverses fractions sont determinees avec le systeme 
10 d'equations suivant : 

0.(x;r ^^)^^Xy^^ -yZ) avec l<i<N 

(i-e-)x-+ery;r = z^ 

et 

N N 

1=1 1=1 

15 

dans lesquelles : 

<Di est la fiigacite du composant i. 

Selon une premiere caracteristique de i'invention, on memorise la 

fraction vaporisee e^et les constantes d'equilibre k"=— ^ avec i variant 

20 entre 1 et N, afin de les utdiser dans une autre etape denommee 
"delimiping" qui sera explicitee ulterieurement. 

Dans une deuxieme etape et pom- chaque maille, on evalue 
certaines des proprietes de chaque phase liquide et gaz, a savoir masse 
volimiique (p™ et p^), viscosite, saturation et permeabilite relative. 

25 Dans une troisieme etape, pour chaque maille j et au pas de 

temps m, on calcule les coefficients de I'equation pour la pression. 

Les equations pour la pression d'une maille j a I'instant t+At ou au 
pas de temps m+1, sont : 
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20 


^Tj 


^pr'-pri_ 


,= ZTTh(pr-pr)+sr 


Pour resoudre ces equations, les coefficients a calculer sont : 
5 - la compressibilite totale C^j (fluide + roche), 

- les transmissivites generalisees 

Trh = Za;r(^oXi+Vi)r'^ih 

M 

10 - le terme source : = J^s™ 

Dans ces equations : 
h sont les mailles voisines de la maiUe j, h e J(j), 

a™ est le volume molaire partiel du composant i dans la maiUe j au pas de 
temps m, X-S et X.^ sont les mobilites des phases huile (o) et gaz (g), 
15 tjh est la transniissivite entre les centres des mailles j et h, 1 etant la maille 

amont. 

Dans une quatrieme etape, on resout pour toutes les mailles les 
eqiiations pour la pression a I'instant t+At c'est-a-dire au pas de temps m+1. 
Ces equations forment un systeme lineaire du type : 


A. p = B 


ou p = (pr*\- - Pr*^ y y J etant le nombre de mailles du reseau. 

Ces equations sont resolues par un precede iteratif standard. 
25 Dans une cinquieme etape, on calcule au pas de temps m+1, les 

debits des phases entre chaque paire de mailles adjacentes et pour chaque 
puits associe a partir, par exemple, de la loi de DARC Y. 

<,=-p-xs.t,H(pr-pr) 

30 
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1 etant la maille amont avec l=j ou h. 

Selon une autre caracteristique de rinvention, on stocke ou 
memorise les debits des phases qui seront utilises dans le delumping. 

Dans une sixieme etape, on evalue la quantite f et la composition 
z du fluide regroupe, dans chaque maille j et au pas de temps m+1 a pardr 
des valeurs au pas de temps m, des debits des phases evalues a la cinquieme 
etape et des compositions des phases calculees a Fetape 1. Pour ce faire, on 
applique les lois de conservation definie par les equations suivantes : 

f;-=f--At|2:k^+"Sh)-sr[ 


r-m-l tn+1 f^y™ — At 


10 Par im processus iteratif, on repete les etapes ci-dessus pour 

calculer les variables principales aux temps t+2At, t-^SAt, etc... ou encore 
aux pas de temps m+1, m+2, m+S, etc., I'increment At pouvant etre 
constant ou variable. 

Selon la presente invention, on defmit une modelisation detaillee 

15 a Q composants et/ou pseudo-composants ainsi qu'une autre equation d'etat 
decrivant le fluide ainsi modelise. 

Pour cela, on choisit prealablement les composants et/ou pseudo- 
composants dont on veut connaitre les caracteristiques et qui sont 
necessaires pour representer le fluide de production dans les conditions de 

20 surface. Le nombre Q qui est superieur au nombre N de la composition 
regroupee est generalement de Tordre de 20 a 30. On connait la composition 
du fluide detaiUe a Q composants et/ou pseudo-composants au pas de temps 
m=0. 

Les donnees stockees aux premiere et cinquieme etapes sont 
25 utilisees pour les etendre aux Q composants et/ou pseudo-composants de la 
composition detaillee. 

Comme I'ecoulement pent s'effectuer dans n'importe quelle 
direction entre les mailles voisines, chaque maille ne peut pas etre traitee 
independamment. C'est pourquoi, a chaque pas de temps, toutes les mailles 
3 0 sont traitees ensemble. 
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Pour realiser le delumping, on fait ITiypothese qu'a chaque pas de 

temps m, 

- la fraction molaire de la phase vaporisee dans chaque maille est 
independante du nombre de constituants du fluide, 

- le nombre de moles de liquide et de vapeur s'ecoulant entre les mailles et 
dans ou hors des puits est independant du nombre de constituants du 
fluide, 

US;h = ur^et USh = u™ 
So7 = s^etSS = s^ 

Les lettres en majuscule sont utilisees pour representer le fluide 
detaille et les lettres en minuscule sont utilisees pour representer le fluide 
regroupe. 

Le calcul des constantes d'equilibre peut etre effectue a partir de 
Texpression : 

Ln[Kj] = AC, +2^^n-^. ou AC„ = Cj; -C^ 

comme indique par C. LEIBOVICI dans son article "A Consistent 
Procedure for Pseudo-Component Delumping" paiu dans FLUID PHASE 
EQUILffiRLA^ 117(1996), 225-232, 

Pour une equation d'etat a deiox parametres, comme celle de 
PENG ROBINSON, Tequation generalisee de LEIBOVICI peut etre ecrite 
pour le fluide detaille de la maniere suivante : 

Ln(K,7) = ii^+P7Va7 + Yrbi 
ou ai et bi sont les parametres de Tequation d'etat pour le composant I et ou 
ies parametres t|",P" ety" sont calcules a partir des constantes d'equilibre 

pour le fluide regroupe (k") en minimisant la fonction : 
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De ce fait, lorsqu'on connait la composition detaillee ZiT dans la 
maille j au pas de temps m, on peut estimer les compositions detaillees des 
phases liquide et vapeur dans la meme maille j et pour le meme pas de 
temps m. 

5 

X™ = r = M T avec 1 < I < Z 

i+0»(k^-i) ' i+0-(k^-i) 

I etant le composant et/ou pseudo-composant de la composition 
detaillee. Ces compositions peuvent etre normalisees afin que leur sonune 
10 soit egale a 1 si necessaire. 

Lorsqu'on connait les compositions des phases detaillees pour 
toutes les mailles du reseau au pas de temps m, il est possible d'estimer les 
compositions detaillees z^' ^e temps consecutif m+1 a partir des 

equations ci-apres : 


15 


ZS'Fr-lYSSS + X^SS)- (Y5?UBh + X5USh) 
■rr'= — avec 1=1,. -.,Q 


OU 


20 Fr'=Fr-(s"+s^)- ^(us.+u5,) 

heJ(j) 

dans lesquelles j -j pour un debit positif de la maille j vers la maille h ou 
dans le puits et j-h pour im debit negatif, j' correspondant au flviide injecte 
pour les puits d'injection ou S sera negatif. ^ 

Grace a la presente invention, on obtient une composition 

25 detaillee a Q composants et/ou pseudo-composants qui permet de mieux 
definir les profils de production previsionnels et ce, avec une precision 
accrue comparee aux methodes anterieures. Cette plus grande precision des 
profils de production previsionnels permet un choix de schema de 
developpement (type et taiUe d'installations de surface, nombre de puits, 

30 etc..) plus fiable, avec pour consequence une augmentation de la rentabilite 
et une diminution du risque economique. 
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Dans ce qui precede, certaines valeurs obtenues au cours de la 
simiilation de gisement sont stockees ou memorisees pour etre utilisees 
ulterieurement dans le delumping. II va de soi que ces memes valeurs 
pourraient etre utilisees directement si on decide d'effecnier le delumping en 
meme temps que ladite simulation de gisement. 
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REVENDICATIONS 


1. Methode de simulation pour predire, en fonction du temps, la 

composition detaillee d'un fluide produit par un reservoir, du type consistant 
a : 

a) representer le reservoir sous la forme d'un reseau de mailles (j) 
dont chacune constitue un volume elementaire rempli de flxiide, 

b) defmir le fluide par une modelisation regroupee a N compo- 
sants et pseudo-composants (i) et a detenniner une equation d'etat decrivant 
le fluide dans cette modelisation regroupee, 

c) definir egalement le fluide par une modelisation detaillee a Q 
composants et/ou pseudo-composants, Q etant superieur a N, et a deter- 
miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa<?on connue en soi, une simulation composi- 
tionnelle du fluide regroupe a N composants et pseudo-composants (i), 
ladite simulation compositionnelle permettant de calculer au mdins pour 
chaque maille (j) et a des pas de temps consecutifs (m, m+1, ...) la fraction 
vaporise (ej*), les constantes d'equilibre liquide-vapeur [k^) de chaque 
composant (i), les debits d'injection ou de production (sj*) et pour chaque 
paire de mailles (j, h) les debits des phases liquide (u^^) et vapeur 

flxiide a N composants et pseudo-composants, et elle est caracterisee en ce 
qu'elle consiste en outre a : 

e) determiner, a chaque pas le temps (m) et pour- chaque maille 
(j) la composition des phases liquide et vapeur du fluide detaille de la 
composition Q du fluide, a partir des valeurs de la fraction vaporisee (e-) et 

des constantes d'equilibre (k^) du fluide regroupe, 

f) evaluer pour chaque maiile, au pas de temps (m+1), la quantite 
molaire de chacxm des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (m), des debits 
des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a I'etape e), et a 
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g) evaluer, pour chaque puits de production, la composition 
detaillee du fluide prodmt entre des instants t et t' correspondant aux pas de 
temps m et m+1, a partir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaiUe determinees a 
5 I'etape e). 

2. Methode selon la revendication 1, caracterisee en ce que les 
etapes e, f et g sont mises en oeuvre en meme temps que la reahsation de la 
simulation compositionnelle du fluide regroupe. 

3, Methode selon la revendication 1, caracterisee en ce que les 
10 resultats de I'etape d sont memorises puis utilises ulteriem-ement pour la 

mise en oeuvre des etapes e, f et g. 
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